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Spectral  Measurements  in  a 
Disturbed  Boundary  Layer  Over  Snow 

EDGAR  L  ANDREAS 


INTRODUCTION 

Atmospheric  turbulence  is  broad-banded.  In  the  surface  layer,  fluctuations  in  the  velocity 
components  and  in  the  scalar  atmospheric  constituents  occur  over  a  continuum  of  frequen¬ 
cies  from  less  than  1/hour  to  well  over  100  Hz.  Computing  spectra  and  cospectra  from  time 
series  of  these  variables  separates  the  contributions  to  the  total  signal  variance  (or  covariance 
in  the  case  of  cospectra)  by  frequency.  Thus,  spectra  and  cospectra  provide  a  picture  of  the 
turbulence  process.  When  topographic  inhomogeneities  perturb  the  atmospheric  boundary 
layer,  such  representations  are  especially  useful  for  examining  how  and  at  what  scale  the  top¬ 
ography  influences  the  turbulence  processes. 

Velocity  and  temperature  spectra  and  cospectra  have  been  measured  often  over  horizon¬ 
tally  homogeneous  land  surfaces  (e.g.,  Kaimal  et  al.  1972,  1976,  Wyngaard  and  Cote  1972, 
Tsvang  et  al.  1985).  Humidity  spectra  have  been  measured  less  frequently  over  land  (Miyake 
and  McBean  1970,  McBean  1971,  Priestley  and  Hill  1985),  because  weak  humidity  fluctua¬ 
tions  make  the  measurements  difficult.  Over  the  ocean,  where  fluctuations  in  humidity  can 
be  large,  spectra  of  the  velocity  components  and  of  the  temperature  and  humidity  have  all 
been  measured  extensively  (e.g..  Pond  et  al.  1971,  Leavitt  1975,  Schmitt  et  al.  1979,  Nicholls 
and  Readings  1981,  Smith  and  Anderson  1984). 

Over  snow-covered  surfaces,  on  the  other  hand,  spectral  measurements  have  been  rare.  In 
fact,  most  spectral  measurements  over  snow  have  been  made  over  snow-covered  sea  ice. 
Banke  and  Smith  (1971,  1973)  and  Andreas  and  Paulson  (1979)  reported  velocity  spectra  and 
cospectra  measured  over  compact  Arctic  sea  ice;  Smith  et  al.  (1970)  reported  similar  meas¬ 
urements  in  the  Gulf  of  St.  Lawrence.  Teunissen  (1980)  measured  spectra  of  the  three  veloc¬ 
ity  components  over  snow-covered  ground  in  farming  country  near  Toronto,  Canada.  Smith 
(1972)  reported  additional  velocity  spectra  from  the  Gulf  of  St.  Lawrence  and  included  a  few 
temperature  spectra  and  cospectra.  But  only  Thorpe  et  al.  (1973)  have  measured  both 
temperature  and  humidity  spectra  simultaneously  over  a  snow-covered  surface — again  sea 
ice.  Hicks  and  Martin  (1972),  MacKay  and  Thurtcll  (1978)  and  Yelagina  et  al.  (1978)  meas¬ 
ured  velocity,  temperature  and  humidity  fluctuations  over  snow-covered  Lake  Mendota  and 
over  snow  fields  in  Ontario,  Canada,  and  near  Leningrad,  respectively,  but  none  reported  in¬ 
dividual  spectra  or  cospectra,  just  covariance  (integrals  of  the  cospectra). 

The  measurements  that  1  report  here  are  thus,  evidently,  the  first  extensive  set  of  turbu¬ 
lence  spectra  and  cospectra  measured  over  snow-covered  ground.  Although  the  experimental 
site  was  fairly  homogeneous  for  a  couple  hundred  meters  around  my  instruments,  hills 
typically  1(X)  m  high  bordered  the  site  on  two  sides.  Hence,  all  the  measured  longitudinal  (u) 
and  vertical  (w)  velocity  spectra  and  the  temperature  (/)  and  humidity  (,q)  spectra  show  iner- 


tial  or  inertial-conveciive  subranges  typical  of  the  horizontally  homogeneous  near-field  but 
contain  excess  energy  at  low  frequency  because  of  far-field  topographic  disturbances.  The 
w-t,  w-q  and  t-q  cospectra  generally  behave  in  the  inertial-convective  subrange  as  cospectra 
collected  over  a  horizontally  homogeneous  surface;  the  «-w  cospectra  usually  have  erratic 
sign  throughout  their  frequency  range,  again  suggesting  just  how  much  the  topography  af¬ 
fects  the  velocity  field.  The  phase  and  coherence  spectra  corroborate  the  topographic  effects 
by  implying  a  dichotomy  between  low-frequency  and  high-frequency  transfer  processes. 

Because  1  have  simultaneous  measurements  of  the  t  and  q  spectra  and  the  t-q  cospectrum, 
it  was  possible  to  make  theoretical  computations  of  refractive  index  spectra  for  light  of  visi¬ 
ble  and  millimeter  wavelengths.  Both  sets  of  refractive  index  spectra  have  inertial-convective 
subranges.  Hence,  lastly,  I  compare  refractive  index  structure  parameters  (Cji)  obtained 
from  the  spectra  for  visible  wavelengths  with  simultaneous  scintillometer  measurements  of 
the  same  quantity. 


MEASUREMENTS 


The  turbulence  data  that  I  will  report  were  all  collected  over  a  field  near  Grayling,  Miciii- 
gan,  that  was  always  covered  with  at  least  30  cm  of  snow.  Figure  1  shows  the  site;  the  wind 
vectors  all  terminate  on  the  instrument  mast.  The  snow  field  was  basically  a  north-south  rec¬ 
tangle  bordered  by  5-m-high  pine  trees.  Between  the  instrument  mast  and  this  tree  line  there 
was  about  900  m  of  field  to  the  north,  300  m  to  the  south,  and  150  m  to  the  west.  In  the 
northeast-to-southeast  sector  were  instrument  vans,  vehicles,  a  few  small  buildings,  and  ac¬ 
cess  roads;  thus,  1  did  not  sample  easterly  winds.  The  hills  in  Figure  1  were  also  forested. 

Figure  2  shows  the  instrument  mast.  The  instruments  for  measuring  the  turbulent  velocity, 
temperature  and  humidity  fluctuations  were  roughly  2  m  above  the  surface.  1  used  TSl  ane¬ 
mometers  and  crossed,  cylindrical  hot-film  sensors  having  lengths  of  1  mm  and  diameters  of 
roughly  53  Mtn  for  measuring  longitudinal  and  vertical  velocity  fluctuations  (Andreas  and 
Murphy  1986b).  Because  these  crossed  sensors  determine  a  vertical  plane,  it  is  important  to 
keep  them  oriented  into  the  wind.  The  instrument  mast,  thus,  could  rotate  around  its  vertical 
axis  so  the  sensors  could  be  pointed  into  the  wind  at  the  start  of  each  recording  run.  I  have 
described  my  procedure  for  converting  the  hot-film  voltage  signals  to  u  and  w  time  series  in 
an  earlier  publication  (Andreas  1979). 

For  measuring  the  temperature  fluctuations  (0,  I  used  a  platinum  resistance  thermometer 
(PRT)  (LaRue  et  al.  1975)  designed  and  built  at  CRREL.  The  sensor,  made  by  DISA  (now 
Dantec  Electronics),  is  a  l-/tm-diameter  platinum  wire  1  mm  long.  For  fluctuations  in  absol¬ 
ute  humidity  (q),  1  used  a  Lyman-alpha  hygrometer  also  built  at  CRREL  and  based  on 
Buck’s  design  (Buck  1976,  1977,  Buck  and  Post  1980).  The  uranium  hydride  source  tube  and 
the  nitric  oxide  detector  tube  in  this  hygrometer  have  a  separation  that  is  variable  from  0  to 
10  cm  for  optimizing  the  response  of  the  hygrometer  to  humidity  fluctuations. 

The  mast  also  held  a  wind  vane  for  orientation  and  instruments  for  obtaining  ambient 
conditions  at  1  and  2  m.  An  R.M.  Young  propeller  vane  yielded  the  wind  speed  (Upv)  and  di¬ 
rection  (0pv)  at  2  m,  and  a  second  propeller  measured  the  wind  speed  at  1  m.  Two  General 
Eastern  systems  measured  average  temperature  (7)  with  a  rugged  PRT  and  and  absolute  hu¬ 
midity  (Q)  with  a  cooled-mirror  dew-point  (frost-point)  hygrometer.  The  T  and  Q  measure¬ 
ments  at  2  m  also  allowed  in  situ  calibration  of  the  Lyman-alpha  hygrometer. 

The  signals  from  the  fast-responding  sensors  were  passed  through  switched-capacitor  fil¬ 
ters  with  cutoff  frequencies  of  1 10  Hz;  the  nominal  mean  of  each  signal  was  then  removed, 
and  the  remaining  fluctuation  amplified  to  fill  a  ±  10-V  range.  All  signals  were  then  digitized 
at  250  Hz  with  12-bit  accuracy  and  written  on  magnetic  tape. 
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Figure  2.  The  instrument  mast,  with  two  levels  for  measuring  the  average  quantities  U,  T  and 
Q,  and  a  2-m  measurement  of  the  turbulent  fluctuations  u,  w,  t  and  q. 

After  the  experiment  1  screened  these  time  series  and  selected  for  further  analysis  the  best 
quality  data — segments  when  the  instruments  were  performing  optimally,  when  conditions 
were  fairly  steady,  and  when  the  wind  direction  provided  the  longest  fetch  over  smooth 
snow.  The  analysis  entailed  removing  linear  trends,  computing  fast  Fourier  transforms 
(FFTs)  of  the  series  in  blocks  of  8192  samples  each,  forming  spectral  densities  from  the 
transform  coefficients,  and  summing  the  spectral  densities  over  logarithmically  spaced  fre¬ 
quency  bands.  In  addition,  I  block-averaged  from  18  to  66  consecutive  8192-sample  blocks 
(depending  on  the  run  length)  to  compute  the  final  band-averaged  spectral  densities.  By  com¬ 
puting  FFTs  of  the  series  of  block  averages,  F  also  extended  the  spectra  to  lower  frequencies 
(Pond  et  al.  1971). 

Because  few  samples  go  into  the  low-frequency  extension  of  the  spectra,  the  spectral  esti¬ 
mates  have  large  uncertainties  here  yet  contain  significant  energy.  To  reduce  the  contribution 
from  these  highly  variable  but  uncertain  spectral  estimates,  I  applied  a  high-pass  digital  RC 
filter  with  a  100-s  time  constant  (Kaimal  and  Haugen  1971)  to  all  the  spectra  and  cospectra. 

During  the  analysis,  I  discovered  that  the  filters  produced  significant  noise  near  their  cut¬ 
off  frequency,  and  this,  naturally,  appeared  in  the  spectra.  Thus,  although  the  raw  spectral 
estimates  should  extend  to  the  Nyquist  frequency,  125  Hz,  because  of  this  filter  noise,  the 
block-  and  band-averaged  spectral  estimates  have  meaningful  information  only  to  about  40 
Hz. 


Table  1.  Summary  of  spectral  runs,  January  1984.  The  average  wind 
speed,  temperature,  humidity  and  wind  direction  were  measured  at  roughly 


2  m. 

Run 

Date 

Time  (EST) 
Start  Stop 

U 

(m  s-') 

T 

rc) 

Q 

(10-'  kg  m  'J 

Direction 

n 

i.i 

10 

1058 

1125 

4.86 

-13.7 

1.07 

1 

1.2 

to 

1125 

1152 

5.33 

-13.3 

0.99 

347 

2.1 

10 

1224 

1255 

4.54 

-13.5 

0.92 

356 

2.2 

10 

1259 

1330 

3.98 

-13.6 

1.10 

5 

3.1 

12 

1027 

1055 

3.00 

-13.1 

0.91 

147 

4.1 

12 

1110 

1146 

4.04 

-14.6 

1.39 

158 

4.2 

12 

1146 

1222 

4.37 

-13.6 

1.40 

162 

5.1 

12 

1245 

1316 

4.56 

-11.6 

1.41 

160 

5.2 

12 

1321 

1352 

4.31 

-11.3 

1.40 

165 

6.1 

14 

1053 

1108 

4.36 

-13.6 

1.28 

30 

6.2 

14 

1109 

1137 

3.92 

-13.5 

1.27 

352 

6.3 

14 

1138 

1206 

3.49 

-13.2 

1.24 

30 

7.1 

14 

1225 

1301 

3.19 

-12.5 

1.26 

25 

7.2 

14 

1301 

1337 

2.40 

-11.9 

1.26 

22 

8.1 

16 

1434 

1448 

3.83 

-  9.2 

1.89 

205 

9.1 

17 

1401 

1428 

3.96 

-  9.5 

1.26 

281 

9.2 

17 

1429 

1456 

4.02 

-  9.3 

1.21 

266 

10.1 

17 

1511 

1526 

3.96 

-  9.5 

1.18 

256 

10.2 

17 

1536 

1548 

4.56 

-  9.7 

1.33 

280 

10.3 

17 

1550 

1621 

4.93 

-  9.6 

1.25 

265 

11.1 

17 

1629 

1702 

4.42 

-11.9 

1.25 

252 

11.2 

17 

1702 

1728 

3.79 

-13.2 

1.29 

246 

12.1 

18 

1318 

1328 

3.62 

-14.2 

1.07 

269 

13.1 

21 

1222 

1236 

4.75 

-17.2 

0.90 

197 

13.2 

21 

1239 

1253 

3.95 

-16.6 

0.91 

215 

14.1 

21 

1328 

1338 

4.27 

-16.7 

0.90 

203 

14.2 

21 

1349 

1411 

4.19 

-16.3 

0.94 

208 

14.3 

21 

1412 

1434 

4.26 

-15.7 

o.% 

203 

15.1 

22 

1423 

1434 

3.62 

-  8.8 

1.99 

208 

Table  1  ii:Ts  the  experi.Tiental  runs  analyzed  and  the  ambient  conditions.  Ail  runs  were 
during  daylight  hours,  and  all  were  from  9  to  36  minutes  in  duration.  The  snow  was  always 
at  least  30  cm  deep;  its  surface  was  usually  smooth  but  occasionally  had  been  sculptured  into 
small  wavelets  by  the  wind.  The  upwind  fetch  over  uniform  snow  depended  on  wind  direc¬ 
tion  but  was  always  150-600  m.  Rarely  was  it  snowing. 

Not  all  runs  in  Table  1  yielded  useful  w,  i  and  q  time  series.  The  delicate  PRT  sensor  broke 
before  or  during  some  runs,  and  the  Lyman-alpha  hygrometer  was  malfunctioning  during 
others.  After  seeing  plots  of  the  spectra  and  cospectra,  I  also  discarded  a  few  because  their 
shapes  were  spurious. 


SPECTRA 


u  and  w  velocity  spectra 

It  is  customary  to  plot  spectra  and  cospectra  nondimcnsionally  so  that  spectra  collected 
under  diverse  conditions  will  collapse  into  universal  form.  In  the  atmospheric  surface  layer, 
that  nondimensionalization  focuses  on  collapsing  velocity  spectra  in  the  inertial  subrange 


and  scalar  spectra  in  the  inertial-convective  subrange.  Here,  with  the  longitudinal  velocity 
spectrum  4>i,  as  an  example,  we  expect  the  Kolmogorov  form: 

*„(*•)  =  a  (1) 

where  k  =  turbulence  wavenumber 

i  =  dissipation  rate  of  turbulent  kinetic  energy 
a  =  Kolmogorov  constant  for  velocity. 

Since  we  measure  frequency  rather  than  wavenumber  spectra,  we  normally  use  Taylor’s  hy¬ 
pothesis,  k  =  2vf/U,  where  /is  the  frequency  and  U  is  the  average  longitudinal  wind  speed 
at  the  measurement  height,  to  convert  eq  I  to  a  nondimensional  frequency  spectrum  (Kaimal 
et  al.  1972): 


/<*>«(/) 


=  (fz/uy 


Thus,  fz/U  is  the  nondimensional  surface  layer  turbulence  frequency,  with  z  the  measure¬ 
ment  height. 

Figure  3  shows  longitudinal  velocity  spectra  made  nondimensional  according  to  eq  2  for 
the  29  runs  listed  in  Table  1 .  Figure  4  shows  the  nondimensional  vertical  velocity  spectra.  As 
1  will  explain  shortly,  e  comes  from  the  spectral  measurements  via  eq  2;  hence  the  inertial- 
subrange  collapse  of  the  spectra  in  Figure  3  is  somewhat  artificial. 

For  comparison  the  dashed  line  in  Figure  3  is  a  slightly  modified  version  of  the  Kansas 
neutral-stability  u  spectrum  (Kaimal  et  al.  1972); 


=  53.9 


[1 +33(/z/t/)]’'’ 


Although  eq  3  has  a  slightly  different  form  than  the  equation  given  by  Kaimal  et  al.  (1972,  eq 
21a),  the  only  real  difference  is  that  eq  3  represents  an  a  value  of  0.54  (Williams  and  Paulson 
1977,  Antonia  et  al.  1981)  instead  of  the  value  of  0.52  that  their  equation  implies.  The  dash¬ 
ed  line  in  Figure  4  is 


=  1.12 


fz/U 

{\+5.l{fz/U)^'^] 


a  modified  version  of  the  Kansas  neutral-stability  w  spectrum  (Kaimal  et  al.  1972,  eq  21c). 
Equation  4  is  different  from  the  form  given  by  Kaimal  et  al.  (1972)  in  that  it  also  reflects  an  a 
value  of  0.54,  and  the  multiplicative  constant  on  the  right  side  is  adjusted  so  that  4>w/4>u  = 
V,  at  large  /,  as  required  by  isotropy.  The  equations  given  by  Kaimal  et  al.  (1972)  yield 
=  ••22  at  large/. 

Neutral-stability  Kansas  spectra  and  cospectra  are  appropriate  for  comparing  with  my 
spectra  and  cospectra  because,  as  1  will  show  shortly,  all  29  runs  in  Table  1  were  made  in 
near-neutral  stability:  the  stability  parameter  z/L,  where  L  is  the  Obukhov  length,  was  al¬ 
ways  between  -0.01  and  0.05.  With  a  measurement  height  of  2  m,  all  spectra  were  collected 
in  the  “inner  region”  defined  by  Jackson  and  Hunt  (1975). 

The  spectra  in  Figures  3  and  4  collapse  well  into  common  shapes,  especially  in  the  inertial 
subrange.  The  scatter  in  the  spectra  at  low  frequency  is  a  consequence  of  both  the  statistical 
uncertainty  in  the  spectral  estimates  at  these  frequencies  and  the  inherent  variability  in  the 
horizontal  wind  at  these  scales. 


loq(f2/U) 


Figure  3.  Nondimensional  u  velocity  spectrum  for  each  run  in  Table  1.  The  dashed  line  is 
the  Kansas  neutral-stability  u  spectrum  (eq  3). 
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Figure  4.  Twenty-two  nondimensional  w  velocity  spectra.  The  dashed  tine  is  the  Kansas 
neutral-stability  w  spectrum  (eq  4). 


These  collected  spectra  are  different  in  two  ways  from  the  Kansas  spectra,  which  were  col¬ 
lected  over  a  horizontally  homogeneous  surface.  They  have  much  more  energy  at  low  fre¬ 
quency,  and  the  inertial  subrange  starts  at  a  nondimensional  frequency  that  is  roughly  a  dec¬ 
ade  lower  than  in  the  Kansas  spectrum.  The  topography  is  responsible  for  both  conditions. 
Whether  the  wind  was  coming  down  from  the  hills  or  being  dammed  as  it  flowed  up  over 
them  (Fig.  1),  the  effect  was  the  same — to  induce  fluctuations  in  the  vertical  velocity.  Be¬ 
cause  of  mass  continuity  and  pressure  correlations,  these  vertical  velocity  fluctuations  creat¬ 
ed  fluctuations  in  the  other  two  velocity  components  at  a  scale  roughly  that  of  the  topog¬ 
raphy.  In  both  Figures  3  and  4  the  maximum  difference  between  the  measured  and  the  Kan¬ 
sas  spectra  occurs  for  log(/z/t/)  between  -2  and  -1.5;  with  z  =  2  m,  this  implies  that  the  ed¬ 
dies  responsible  for  this  energy  had  wavelengths  of  60-l(X)  m,  about  the  height  of  the  hills 
above  the  snow  field.  It  is  obvious  from  Figure  1  that  there  are  a  host  of  other  longer  and 
shorter  length  scales  associated  with  the  topography,  but  the  height  of  the  hills  is  the  simplest 
and  the  one  most  clearly  represented  in  the  velocity  spectra  (e.g.,  Jackson  and  Hunt  1975). 

The  second  difference  in  my  spectra,  the  extended  inertial  subrange,  is  likewise  a  conse¬ 
quence  of  the  topography.  Because  of  the  topographic  energy  generation,  the  spectral  peaks 
in  Figures  3  and  4  appear  at  lower  frequencies  than  in  the  Kansas  spectra.  In  other  words,  be¬ 
cause  of  the  topography,  the  energy-containing  or  energy-producing  region  is  shifted  to 
longer  length  scales,  thereby  increasing  the  spectral  separation  between  the  production  and 
dissipation  regions.  The  inertial  subrange  can,  therefore,  begin  at  a  proportionately  lower 
nondimensional  frequency  (Tennekes  and  Lumley  1972,  p.  264ff.). 

Although  a  few  measurements  of  velocity  spectra  exist  for  complex  terrain,  the  terrain  at 
this  site  was  more  complex  than  most.  Not  only  was  there  a  series  of  hills,  rather  than  the  iso¬ 
lated  hill  commonly  studied  (Jackson  and  Hunt  1975,  Bradley  1980,  Smedman  and  Berg¬ 
strom  1984,  Mason  and  King  1985),  the  surface  roughness  (zo)  of  the  site  was  small  (~0.1-I 
cm)  (Banke  et  al.  1980)  over  the  near-field  snow  but  relatively  large  ( -0.2-0.5  m)  (Garratt 
1977)  over  the  wooded  hills  and  the  woods  surrounding  the  snow  field.  My  velocity  spectra, 
nevertheless,  show  some  of  the  same  trends  seen  by  others. 

For  example,  though  Smedman  and  Bergstrom  (1984)  made  measurements  at  the  summit 
of  a  fairly  low,  isolated  hill,  their  velocity  spectra  from  the  inner  region  collected  under  near¬ 
neutral  stability  had  characteristics  like  mine.  Both  their  u  and  w  spectra  behaved  as  at 
high  frequency,  with  <!>*  somewhat  lower  than  the  Kansas  neutral-stability  spectrum  here, 
like  my  w  spectra  (Fig.  4).  At  low  frequency,  their  u  and  w  spectra  contained  more  energy 
than  the  Kansas  spectra  though  not  as  much  excess  energy  as  mine  do,  perhaps  because  their 
hill  was  only  20  m  above  the  surrounding  terrain. 

Mason  and  King  (1985)  measured  u  and  w  spectra  in  neutral  stability  conditions  within  the 
inner  region  at  the  summit  and  near  the  base  of  an  isolated,  100-m-high  hill — a  height  typical 
of  the  topography  for  this  study.  They,  too,  found  excess  energy  at  low  frequency  compared 
to  the  Kansas  spectra  for  both  velocity  components  and  at  both  locations.  At  high  fre¬ 
quency,  however,  both  u  and  w  spectra  rolled  off  faster  than  In  an  earlier  study  of  flow 
over  ridges  and  valleys  using  basically  the  same  instruments.  Mason  and  King  (1984),  how¬ 
ever,  did  find  /  behavior  in  both  u  and  w  spectra  at  a  ridge  crest  and  on  the  valley  floor. 
And  again  they  found  excess  energy  in  <!>„  and  4>w  at  low  frequency,  with  spectra  from  the 
valley  floor  containing  more  energy  than  spectra  from  the  ridge  crest. 

Teunissen  (1980)  measured  spectra  of  the  three  velocity  components  under  near-neutral 
conditions  over  relatively  flat,  snow-covered  terrain  studded  with  occasional  trees,  fences 
and  farm  buildings.  Like  mine,  his  spectra  showed  more  energy  than  the  Kansas  neutral- 
stability  velocity  spectra;  he  attributed  this  to  the  rough  terrain  at  his  site.  Teunissen  also  no¬ 
ticed  that  the  low-frequency  behavior  of  his  spectra  made  them  look  more  like  Kansas  spec¬ 
tra  collected  in  unstable  conditions.  He,  however,  discounted  the  notion  that  unstable  strati¬ 
fication  higher  in  the  boundary  layer  affected  his  spectra,  because  during  his  experiment  the 
boundary  layer  had  slightly  stable  stratification.  He,  thus,  reiterated  his  conclusion  that  the 


roughness  of  his  site,  not  the  stratification,  led  to  the  excess  low-frequency  energy.  I  can 
make  a  similar  argument.  Radiosondes  released  during  our  experiment  typically  showed  an 
inversion  layer  to  no  higher  than  300  m,  with  a  well-mixed  layer  usually  above  this.  Thus,  al¬ 
though  the  surface  layer  was  nearly  neutral,  conditions  higher  up  were  slightly  stable;  there¬ 
fore,  the  excess  energy  I  find  at  low  frequency  in  the  velocity  spectra  cannot  be  a  conse¬ 
quence  of  unstable  stratification  above  the  surface  layer. 

At  the  Boulder  Atmospheric  Observatory  (BAO),  a  site  with  much  gentler  local  relief  than 
mine,  Kaimal  et  al.  (1982)  measured  velocity  spectra  in  unstable  conditions  that  were  much 
the  same  as  the  Kansas  spectra.  The  only  differences  they  noted  were  that  the  BAO  w  spectra 
had  peaks  shifted  to  somewhat  lower  frequency,  and  consequently  these  spectra  contained 
slightly  more  energy  at  low  frequency  than  w  spectra  collected  over  horizontally  homogen¬ 
eous  terrain.  At  BAO  in  stable  conditions,  on  the  other  hand.  Hunt  et  al.  (1985)  saw  the  ef¬ 
fects  of  topographically  induced  internal  gravity  waves  at  low  frequency  in  w  spectra  collect¬ 
ed  near  the  surface.  At  high  frequency  the  inertial  subrange  seemed  unaffected. 

In  summary,  my  measured  u  and  w  spectra  are  consistent  with  most  other  velocity  spectra 
measured  in  the  inner  region  over  complex  terrain.  Both  u  and  w  spectra  have  extensive  iner¬ 
tial  subranges  where  the  small  scales  have  come  rapidly  to  equilibrium  with  the  underlying 
surface  and  the  turbulence  seems  locally  isotropic  (Panofsky  et  al.  1982,  Mason  and  King 
1985).  At  lower  frequency,  in  near-neutral  stability,  the  topography  is  responsible  for  in¬ 
creasing  the  energy  content  in  both  velocity  components  compared  to  spectra  obtained  over 
horizontally  homogeneous  surfaces. 

Temperature  and  humidity  spectra 

At  large  Reynolds  numbers  an  inertial-convective  subrange  should  exist  in  scalar  spectra 
where,  with  the  temperature  spectrum  4>t  as  an  example  (Corrsin  1951,  Phelps  and  Pond 
1971,  Tennekes  and  Lumley  1972), 

4-, (it)  =  (5) 

Here  0i  is  the  Kolmogorov  constant  for  temperature,  and  Nf  is  the  dissipation  rate  of  half  the 
temperature  variance.  Because  the  velocity  field  advects  scalar  contaminants,  however,  we 
would  expect  scalar  spectra  to  mirror  the  energy  levels  of  velocity  spectra  at  low  frequency. 
Despite  low  temperatures,  a  frozen  surface  and  near-neutrality  in  the  surface  layer.  Figures  5 
and  6  show  that  the  turbulent  temperature  and  humidity  fields  during  the  experiment  con¬ 
tained  enough  energy  to  yield  excellent  spectra.  And  both  spectral  sets  have  the  expected 
features — an  inertial-convective  subrange  and  high  energy  at  low  frequency. 

The  dashed  line  in  both  figures  is  a  modified  version  of  the  Kansas  neutral-stability  tem¬ 
perature  spectrum  (Kaimal  et  al.  1972), 

=  34.8  - -  ,  fz/U  <0.15,  (6a) 

(1 -t-24(/z/(/)]’^’ 

=  15.9  - -  ,  fz/U  2  0.15.  (6b) 

[l  +  12.5(/-z/t/)l’'’ 

Although  this  form  looks  different  from  eq  2Id  in  Kaimal  et  al.  (1972),  the  only  essential  dif¬ 
ferences  are  that  eq  6  reflects  the  more  commonly  accepted  ratio  for  the  turbulent  diffusivi- 
ties  for  heat  (Kh)  and  momentum  (Km)  at  neutral  stability  of  Kh/K^  =  1  rather  than  the  val¬ 
ue  0.74  used  by  them,  and  it  implies  |(i,  =  0.80  instead  of  the  value  0.82  implicit  in  their  equa- 
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Figure  5.  Twenty-four  nondimensional  temperature  spectra.  The  clashed  line  is  the  Kansas 
neutral-stability  temperature  spectrum  (eq  6). 
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Figure  6.  Twenty-four  nondimensional  humidity  spectra.  The  dashed  line  is  the 
Kansas  neutral-stability  temperature  spectrum  (eq  6). 


The  dashed  line  in  Figure  6  is  also  given  by  eq  6,  since  Kaimal  et  al.  (1972)  measured  no  hu¬ 
midity  spectra.  In  Figure  6  the  absolute  humidity  spectrum,  4>q,  is  made  nondimensional 
with  A/q,  the  dissipation  rate  of  half  the  humidity  variance.  Notice  that  in  using  eq  6  to  also 
represent  humidity,  I  am  assuming  that  the  Kolmogorov  constant  for  humidity,  /jq,  is  equal 
to  (3,,  =  (3q  =  ^  =  0.80,  as  many  measurements  have  suggested  (Paquin  and  Pond  1971, 

Wyngaard  and  Cote  1971,  Champagne  et  al.  1977,  Raupach  1978,  Wyngaard  et  al.  1978, 
Bradley  et  al.  1981,  Large  and  Pond  1982). 

As  anticipated,  at  high  frequency  the  collected  temperature  and  humidity  spectra  in  Fig¬ 
ures  5  and  6  have  the  behavior  indicative  of  an  inertial-convective  subrange.  At  low  fre¬ 
quency,  although  stability  conditions  were  nearly  neutral,  the  scalar  spectra  contain  roughly 
an  order  of  magnitude  more  energy  than  predicted  by  the  Kansas  measurements.  This  excess 
energy  is  presumably  again  a  consequence  of  the  topography;  large-scale,  topographically  in¬ 
duced  velocity  fluctuations  simply  advect  the  imbedded  scalar  field,  giving  it  spectral  charac¬ 
teristics  similar  to  the  longitudinal  velocity  spectrum.  With  this  excess  spectral  energy  at  low 
frequency,  the  energy-producing  or  energy-containing  region  of  the  scalar  spectra  is  shifted 
to  lower  frequencies,  as  noted  with  the  velocity  spectra.  Hence,  the  inertial-convective 
subrange  seems  to  start  at  lower  frequency  in  the  measured  t  and  q  spectra  than  in  the  Kansas 
spectrum. 

Because  these  temperature  and  humidity  spectra  have  similar  shapes,  it  is  reasonable  to 
compare  both  to  the  Kansas  temperature  spectrum  (Schmitt  et  al.  1979).  The  only  obvious 
difference  between  the  nondimensional  spectra  is  that  the  /  spectra  generally  contain  slightly 
more  energy  at  low  frequency  than  the  q  spectra.  Consequently  the  peaks  in  the  t  spectra  are 
near  \og(Jz/U)  =  -2.3,  while  for  the  q  spectra  they  are  roughly  at  log(/z/t/)  =  -  1.7. 

Rarely  have  simultaneous  measurements  of  temperature  and  humidity  fluctuations  yielded 
such  compatible  spectral  sets.  Most  simultaneous  t  and  q  measurements  have  been  made  over 
the  ocean.  Phelps  and  Pond  (1971)  found  significant  differences  between  temperature  and 
humidity  spectra  measured  during  BOMEX  (Barbados  Oceanographic  and  Meteorological 
Experiment)  and  again  over  the  ocean  near  San  Diego.  In  particular,  their  BOMEX  t  spectra 
had  no  inertial-convective  subrange.  Nicholls  and  Readings  (1981)  also  observed  different 
spectral  forms  for  temperature  and  humidity  obtained  from  an  aircraft  flying  in  the  surface 
layer  above  the  ocean.  Their  humidity  spectra  were  generally  compatible  with  the  Kansas 
temperature  spectra;  but,  although  their  t  spectra  showed  an  inertial-convective  subrange, 
they  had  local  minima  where  the  neutral-stability  Kansas  spectrum  had  its  peak.  Leavitt’s 
(1975)  t  and  q  spectra  from  BOMEX  had  similar  inertial-convective  subrange  behavior,  but 
his  q  spectra  generally  had  much  more  energy  at  low  frequency  than  his  t  spectra.  Friehe  et 
al.  (1975)  found  similar  t  and  q  spectra  over  the  freshwater  Salton  Sea;  but  over  the  ocean 
southwest  of  San  Diego,  their  measured  q  spectra  had  an  extensive  inertial-convective  sub¬ 
range,  while  their  t  did  not.  Schmitt  et  al.  (1979)  measured  humidity  spectra  in  the  North  Pa¬ 
cific  that  generally  looked  like  the  Kansas  temperature  spectra  except  for  an  unexplained 
hump  near/j/t/  =  1  and  excess  energy  at  low  frequency  in  stable  conditions,  which  they  at¬ 
tributed  to  instrument  problems.  Their  temperature  spectra,  however,  were  anomalous. 
These  latter  measurements  and  those  by  Friehe  et  al.  (1975)  help  explain  some  of  the  more 
conspicuous  differences  between  observed  over-ocean  temperature  and  humidity  spectra. 
Sea  salt  particles  can  contaminate  small,  fast-responding  temperature  sensors  by  making 
them  also  responsive  to  humidity  fluctuations  (Schmitt  et  a).  1978,  1979,  Fairall  et  al.  1979). 
Spurious  temperature  spectra  often  result. 

No  similar  complications  exist,  of  course,  with  over-land  turbulence  measurements.  Over 
land,  however,  humidity  fluctuations  may  be  so  small  that  only  sensors  with  very  low  noise 
can  resolve  q  spectra.  For  example,  McBean’s  (1971)  temperature  spectra  showed  a  distinct 
inertial-convective  subrange;  but  his  humidity  spectra,  though  having  levels  similar  to  his  i 


m 
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II 


spectra  at  low  frequency,  fell  off  much  more  rapidly  than  at  high  frequency.  Kohsiek 
(1984)  made  similar  observations;  his  t  spectra  went  as  at  high  frequency,  but  his  q 
spectra  fell  off  somewhat  more  rapidly.  The  temperature  and  humidity  spectra  Schmitt  et  al. 
(1979)  measured  over  Minnesota  farmland  behaved,  generally,  like  the  Kansas  temperature 
spectra;  but,  as  with  their  oceanic  measurements,  the  humidity  spectra  had  an  unusual  hump 
near/z/(y  =  1 .  Although  the  spectra  measured  by  Priestley  and  Hill  (1985)  did  not  extend  to 
frequencies  below  the  spectral  peak,  their  temperature  and  humidity  spectra,  like  mine, 
showed  a  distinct  inertial-convective  subrange  that  seemed  to  begin  at  frequencies  lower 
than  suggested  by  the  Kansas  spectra.  The  t  and  q  spectra  that  Ohtaki  (1985)  measured  over 
vegetated  fields  in  Japan  also  had  an  inertial-convective  subrange;  but  his  t  spectra  had 
somewhat  more  low-frequency  energy  than  his  q  spectra,  as  with  my  spectra. 

The  only  other  temperature  and  humidity  spectra  measured  over  snow  that  1  know  of, 
those  by  Thorpe  et  al.  (1973)  from  the  Beaufort  Sea,  are  dissimilar.  Their  humidity  spectra 
generally  had  an  inertial-convective  subrange,  but  their  temperature  spectra  did  not  because 
of  an  inadequacy  in  their  temperature  sensor.  At  low  frequency  their  temperature  spectra 
had  slightly  more  energy  than  their  humidity  spectra,  as  mine  do. 

Very  few  measurements  of  scalar  spectra  have  been  made  over  complex  terrain.  At  the 
Boulder  Atmospheric  Observatory,  where  the  terrain  is  gently  rolling,  Kaimal  et  al.  (1982) 
found  no  essential  differences  between  measured  temperature  spectra  and  the  Kansas  spectra 
for  unstable  conditions.  In  stable  conditions  at  the  same  site,  however.  Hunt  et  al.  (1985)  saw 
evidence  at  low  frequency  in  their  temperature  spectra  of  internal  gravity  waves  that  they  im¬ 
ply  were  topographically  induced. 

To  summarize,  my  measured  temperature  and  humidity  spectra  are  comparable  with  other 
measured  scalar  spectra.  Their  main  distinguishing  feature  is  the  excess  energy  at  low  fre¬ 
quency,  a  consequence  of  the  topographically  induced  velocity  fluctuations  at  these  scales. 
Because  of  this  excess  energy,  the  spectral  production  region  is  at  lower  frequency  than  in 
spectra  collected  over  a  horizontally  homogeneous  surface;  hence,  the  inertial-convective 
subrange  begins  at  lower  frequency.  This,  too,  has  been  observed  before.  That  t  spectra  have 
somewhat  more  low-frequency  energy  than  q  spectra  is  not  a  hard-and-fast  rule  but  also  has 
been  seen  before. 


INERTIAL-DISSIPATION  ESTIMATES 

With  spectral  measurements  of  velocity  and  scalar  fluctuations  that  show  inertial  or  iner¬ 
tial-convective  subranges,  it  is  reasonable  to  compute  inertial-dissipation  estimates  of  the 
turbulent  momentum  (r)  and  sensible  (//j)  and  latent  (//, )  heat  fluxes: 


T  =  -  QUW  =  Q  Ui 

(7) 

=  QCpWl  =  -QCpU.I, 

(8) 

=  L^wq  =  -L^u.q, 

(9) 

where  q  =  air  density 

Cp  =  specific  heat  of  air  at  constant  pressure 
Z-s  =  latent  heat  of  sublimation  of  ice 
«„  t,  and  q,  -  velocity,  temperature  and  humidity  scales  defined  as  shown 
and  an  overbar  indicates  a  time  average. 
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The  inertial-dissipation  technique  starts  with  budget  equations  for  the  turbulent  kinetic 
energy  and  for  the  scalar  variances.  Fairall  and  Larsen  (1986)  discussed  these  equations  and 
explained  simplifications  that  make  them  tractable.  For  kinetic  energy  the  budget  equation 
seems  to  involve  different  terms  depending  on  whether  conditions  are  stable  or  unstable;  but, 
for  conditions  near  neutral,  all  of  the  equations  reduce  essentially  to  a  balance  between  pro¬ 
duction  and  dissipation.  Specifically, 

{,u\/kz)[<t>m(z/L)~z/L]  =  (  (10) 

[{i^ir/u.Hz]Oy,{z/L)  =  N^  (11) 

((vv^)Vm.i<j]  OhIj/Z-q)  =  A/q  (12) 

where  x  =  von  Karman’s  constant  (0.4) 

=  dimensionless  longitudinal  velocity  gradient 
vj,  =  dimensionless  scalar  gradient. 

Large  and  Pond  (1981,  1982)  used  these  same  equations  for  all  stabilities.  It  is  customary  to 
assume  that  Oh  has  the  same  form  for  both  potential  temperature  and  absolute  humidity 
(Fairall  and  Larsen  1986),  since  this  gradient  function  has  rarely  been  measured  for  humidity 
(Dyer  1974,  Dyer  and  Bradley  1982).  Notice,  however,  in  eq  11  and  12,  1  use  different  stabili¬ 
ty  parameters  for  temperature  and  humidity: 


L 

(13) 

^  Ui  c-)-o.6igl 

(14) 

where  L  is  the  usual  Obukhov  length  corrected  for  moisture  flux,  and  Lq  is  a  length  scale 
based  on  the  moisture  flux  alone  (Busch  1973).  f  andQ  are  representative  surface  layer  val¬ 
ues  of  potential  temperature  and  absolute  humidity.  Later,  we  will  see  that  the  temperature 
and  humidity  fluctuations  were  seldom  positively  correlated  for  the  measurements  reported 
here.  McBean  (1971)  explained  that  unless  a  passive  scalar  contaminant  is  well  correlated 
with  temperature,  L  is  not  the  proper  scaling  length  for  it.  Rather,  a  length  scale  based  only 
on  the  flux  of  the  scalar,  such  as  Lq  is  for  humidity,  is  a  more  meaningful  parameter  (An¬ 
dreas  and  Makshtas  1985). 

A  plethora  of  forms  for  Om  and  exists  in  the  literature  (e.g..  Dyer  1974,  Yamamoto 
1975,  Yaglom  1977).  The  representations  1  use  in  eq  10-12  are  the  same  ones  that  Large  and 
Pond  (1982)  used; 


for  unstable  conditions  (z/L,  z/Lq  <  0), 

o^(z/L)  =  [1-16(t/L)1  '  “  (15) 

Oh(z/Z.)  =  [l-16(c/L)]  '  -  (16) 


-1/2 


Ohl^/Z-g)  -  [1-16(Z/Z.q)1 


(17) 


for  stable  conditions  (z/L,  z/Lq  >  0), 
<t>^(z/L)  =  H(z/L)  =  1  +  7(z/Z.) 


(18) 


<t>iAz/LQ)  =  \+1(z/Lq).  (19) 

The  method  of  estimating  t,  N^  and  for  use  in  eq  10-12  defines  the  inertial-dissipation 
method.  I  found  e  from  measurements  of  <t>i,  in  the  inertial  subrange  using  eq  1  with  a  - 
0.54.  Similarly,  knowing  t  now,  I  found  Ni  and  Nq  from  measurements  of  4>t  and  <I>q  in  the 
inertial-convective  subrange  using  eq  5  for  temperature, 

4>q(A:)  = /iq7Vqr'^3^-5'3  (20) 

for  humidity,  and  with  /3,  =  |3q  =  0.80.  Having  these  values,  I  iteratively  solved  eq  10-12  for 
u„  ^  and  wq. 

Notice  that  eq  1 1  and  12  cannot  yield  the  signs  of  wt  and  wq.  To  determine  these  1  looked 
at  the  inertial-subrange  behavior  of  the  w-t  and  w-t  cospectra  and  at  the  difference  between 
temperatures  and  humidities  measured  at  the  snow  surface  and  at  2  m.  The  surface  tempera¬ 
ture  and  humidity  measurements  were  based  on  a  newly  developed  hygrometric  technique 
(Andreas  1986). 

Although  all  of  the  computations  involved  in  the  inertial-dissipation  technique  are  formal¬ 
ly  possible  because  <l>u,  4>,  and  <I>q  have  /  regions,  I  must  nevertheless  discuss  their  theor¬ 
etical  justification.  The  assumptions  that  went  into  eq  10-12  are  that  the  surface  is  horizon¬ 
tally  homogeneous — thus,  advection  is  unimportant — and  that  turbulent  transport  is  negligi¬ 
ble  (Fairall  and  Larsen  1986).  Even  if  the  turbulent  transport,  which  has  never  been  investi¬ 
gated  over  complex  terrain,  is  negligible,  we  have  already  seen  how  advection  of  the  scalar 
fields  by  the  topographically  disturbed  velocity  field  affects  the  spectra.  But  these  effects 
were  confined  to  long  wavelengths— wavelengths  typical  of  the  topography.  Thus,  as  Panof- 
sky  et  al.  (1982)  and  Mason  and  King  (1985)  suggested,  over  complex  terrain  the  small-scale 
turbulence  still  rapidly  reaches  equilibrium  with  the  local  surface.  Hence,  eq  1,  5  and  20 
should  yield  fairly  accurate  estimates  of  f,  N,  and  Nq.  Although  horizontal  homogeneity  is 
also  the  basis  for  the  flux-gradient  representation  of  the  vertical  profiles — (i>^(z/L),  <^>h(^/Z.) 
and  (bhiz/Lq)  in  eq  10-12— Smith  et  al.  (1983)  and  Andreas  and  Murphy  (1986a)  have  recent¬ 
ly  shown  that  the  functions  (eq  15-19)  are  accurate  and  useful  even  over  extremely  inhomo¬ 
geneous  surfaces.  Consequently,  in  light  of  our  need  for  estimates  of  the  turbulent  fluxes— 
but  keeping  the  possible  shortcomings  of  eq  10-12  in  mind — it  seems  justifiable  to  base 
rough  estimates  of  the  fluxes  on  eq  10-12. 

Table  2  lists  the  results  of  the  inertiai-dissipation  calculations.  A  positive  flux  is  upward — 
from  the  snow  into  the  air.  Where  there  are  blanks  in  the  table,  the  temperature  or  humidity 
spectrum  was  unavailable;  to  make  the  calculations  for  these  runs,  I  assumed  that  the  cor¬ 
responding  dissipation  rate  was  zero.  This  assumption  had  little  effect  on  the  resulting  values 
since  the  heat  fluxes  were  so  small  and  the  absolute  values  of  L  and  Lq  were  so  large. 

Table  2  shows  that  stability  conditions  were  nearly  neutral:  z/L  was  always  between  -0.01 
and  0.05.  The  heat  fluxes  and  were  small,  with  absolute  values  seldom  above  10  W 
m  ^  This  is  not  surprising,  considering  the  gentle  winds  (Table  1)  and  the  heat  fluxes  over 
snow  and  snow-covered  sea  ice  that  others  have  measured  (Hicks  and  Martin  1972,  Thorpe  et 
al.  1973,  McKay  and  Thurtell  1978,  Yelagina  et  al.  1978,  Andreas  and  Makshtas  1985).  One 
interesting  result  evident  in  Table  2  is  that  and  Hy  generally  had  opposite  signs.  Even  dur¬ 
ing  daylight  hours  //,  was  usually  downward,  while  Hy  was  predominantly  upward.  That  is, 
the  air  was  usually  losing  sensible  heat  to  the  surface  but  was  recouping  some  of  this  loss  as 
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Table  2.  Inertial-dissipatioii  estimates  of  the  fluxes. 

L  L>q  u« 


Run 

z/L 

(m) 

(m) 

(m  $-') 

(Wm") 

fIV  m  V 

1.1 

0.011 

190 

-  9800 

0.27 

-  9.5 

3.1 

1.2 

0.008 

250 

-  8900 

0.28 

-  7.8 

3.8 

2.1 

0.011 

180 

-  6400 

0.24 

-  7.0 

3.3 

2.2 

0.017 

120 

-  4000 

0.20 

-  6.5 

3.4 

3.1 

-0.005 

-390 

-11000 

0.19 

1.6 

1.0 

4.1 

-0.003 

-600 

-  5000 

0.21 

1.2 

3.0 

4.2 

0.003 

610 

-  4600 

0.21 

-  1.5 

3.2 

5.1 

0.003 

680 

-  3400 

0.24 

-  2.1 

6.1 

5.2 

0.003 

750 

-  3000 

0.23 

-  1.7 

5.9 

6.1 

0.017 

120 

0.19 

-  5.4 

6.2 

0.022 

92 

0.18 

-  5.9 

6.3 

0.022 

91 

0.17 

-  4.9 

7.1 

0.027 

72 

0.14 

-  3.6 

7.2 

0.050 

40 

0.11 

-  3.2 

8.1 

-  5200 

0.23 

3.5 

9.1 

0.004 

490 

-  4300 

0.22 

-  2.3 

4.1 

9.2 

0.005 

360 

-  3800 

0.21 

-  2.6 

3.9 

10.1 

0.006 

330 

-  3600 

0.22 

-  2.9 

4.4 

10.2 

0.003 

610 

-  6600 

0.25 

-  2.5 

3.7 

10.3 

0.007 

270 

-  5000 

0.26 

-  5.9 

5.3 

11.1 

0.015 

130 

-19000 

0.22 

-  7.7 

0.9 

11.2 

0.026 

74 

6400 

0.20 

-  8.9 

-1.9 

12.1 

-0.005 

-430 

-  3400 

0.21 

1.6 

4.1 

13.1 

-0.008 

-240 

-33000 

0.27 

7.5 

1.0 

13.2 

0.013 

150 

-24000 

0.26 

-10.2 

1.2 

14.1 

-0.010 

-190 

-15000 

0.26 

8.1 

1.8 

14.2 

-16000 

0.26 

1.7 

14.3 

-16000 

0.27 

1.9 

15.1 

-  5000 

0.22 

3.3 

latent  heat  from  surface  sublimation.  Most  other  measurements  of  //l  over  a  snow-covered 
surface  have  also  shown  the  predominance  of  sublimation  over  condensation  (McKay  and 
Thurtell  1978,  Yelagina  et  al.  1978,  Andreas  and  Makshtas  1985). 

The  next  section  discusses  measurements  of  the  w-t  and  w-q  cospectra.  Integrals  of  these 
cospectra  yield  eddy-correlation  estimates  of  //,  and  //l.  Although  these  eddy-correlation  es¬ 
timates  are  broadly  compatible  with  the  inertial-dissipation  estimates  of  and  Wl,  they 
have  much  more  scatter  because  of  large  but  statistically  uncertain  contributions  to  the  inte¬ 
grals  from  the  low-frequency  end  of  the  cospectra.  Hence,  the  inertial-dissipation  estimates 
of  the  fluxes  are  more  useful,  especially  for  nondimensionalizing  the  cospectra. 


COSPECTRA 

With  simultaneous  u,  w,  t  and  q  time  series,  it  is  possible  to  compute  cospectra  (Co^y)  and 
quadrature  spectra  (Qd^y)  of  any  two  variables  Jr  and  y.  The  cospectra,  individually,  are  in¬ 
teresting;  but  since  the  quadrature  spectra  contain  basically  phase  information,  they  are 
most  useful  for  defining  phase  and  coherence  spectra.  If  and  <l>y  are  spectra  of  jt  and  y 
fluctuations,  the  phase  and  coherence  spectra  are  defined,  respectively,  as 


OxyC/)  =  arctan(Q<f,y(/")/Co^y(/')] 


(21) 


(22) 


Ch^yif)  = 


co,y(n’  +  grfxycn^ 


1>x(/')4-y(/^) 


The  phase  spectrum  gives  the  phase  lag  between  simultaneous  fluctuations  in  x  and  y.  The 
coherence  spectrum  is  like  a  spectral  correlation  coefficient:  high  values  of  Ch^y  mean  at  and 
y  have  good  positive  or  negative  correlation;  low  values  mean  they  are  poorly  correlated. 
Through  most  of  the  spectrum,  C/i^y  <  1.0(Panofsky  and  Dutton  1984,  p.  73).  But  because 
of  the  averaging  used  in  the  extended  portion  of  my  spectra  and  cospectra,  Ch^y  can  be  great¬ 
er  than  one  here  as  a  consequence  of  the  statistical  variability  in  these  low-frequency  esti- 


Notice  that  the  phase  and  coherence  spectra  allow  alternative  definitions  of  the  cospec¬ 
trum  and  the  quadrature  spectrum: 


COxy(/)  =  Gxy(/)|  COS[exy(/)l 


Qd„y{f)  =  jGxy(/-)i  sin[exy(/)] 


Gxy(/)l  =  (<l>x(/')<i>y(/^)C/lxy(/)l' 


Thus,  the  quadrant  in  which  the  phase  spectrum  lies  is  directly  related  to  the  signs  of  the  co¬ 
spectrum  and  quadrature  spectrum. 


u-w  cospectra 

The  u-w  cospectral  plots  show  best  the  effects  of  the  topography  on  the  surface  boundary 
layer  at  our  site.  Figure  7a  contains  u-w  cospectra  plotted  nondimensionally  with  typical 
surface-layer  scaling  (Kaimal  et  al.  1972)  for  the  29  runs  in  Table  1.  The  vacillating  sign  of 
the  cospectral  estimates  at  both  high  and  low  frequency  is  obvious.  Over  a  horizontally  ho¬ 
mogeneous  surface,  Cou*  would  be  uniformly  negative — the  wind  transferring  momentum 
downward  to  the  surface  at  all  scales.  The  variability  in  Figure  7a  implies,  however,  that  at 
some  scales  the  air  was  transferring  momentum  downward,  while  at  other  scales  the  transfer 
was  upward.  The  u-w  phase  and  coherence  spectra  (Fig.  7b  and  c)  show  similarly  confused 
pictures.  Oy*  and  C/i,,*  can  have  virtually  any  values  at  low  frequency.  That  is,  sometimes  u 
and  w  were  fairly  well  correlated  here,  but  at  other  times  they  were  poorly  correlated.  At 
high  frequency,  ©uw  >s  generally  between  0°  and  180°  and  decreases  to  zero.  There¬ 
fore,  in  the  inertial  subrange  the  u  fluctuations  often  led  the  w  fluctuations  but  seldom  by 
enough  to  make  COu*  negative. 

Others  have  observed  sign  reversals  in  u-w  cospectra  or  in  time  series  of  the  product  uw; 
but  rarely  have  these  been  as  prevalent  as  in  Figure  7a.  Smith  et  al.  (1970),  for  example, 
found  a  few  positive  cospectral  values  at  the  lowest  frequencies  in  their  measurements  near 
the  surface  of  broken  sea  ice  in  the  Gulf  of  St.  Lawrence.  Pond  et  al.  (1971)  observed  occa¬ 
sional  positive  values  of  the  ocean  for  fz/U  <  0.01  and  attributed  these  to  statisti¬ 

cal  variations.  McBean  and  Miyake  (1972)  saw  positive  values  of  Co^w  at  still  higher  fre¬ 
quencies  in  stable  conditions  over  a  site  that  was  horizontally  homogeneous  on  a  microscale, 
at  least,  and  inferred  that  internal  waves  caused  this  behavior.  Finally,  Haugen  et  al.  (1971) 
reported  occasions  during  the  Kansas  experiment  when  the  uw  product  series  was  positive  as 
a  consequence,  probably,  of  horizontal  roll  vortices  (Brown  1980). 
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c.  Coherence  spectra. 


Figure  7  (cont’d).  Plots  of  the  u-w  cospectrum,  phase  spectrum  and  coherence 
spectrum  for  each  of  the  runs  in  Table  I. 


The  influence  of  the  topography  on  the  u-w  cospectra  becomes  more  obvious  when  we 
compare  Figures  7  and  8.  Figure  8a  shows  six  m-w  cospectra  that  are  negative  throughout  the 
inertial  subrange.  The  runs  that  produced  these  cospectra,  2.2,  6.1,  6.2,  6.3,  7.1  and  7.2 
(Table  1),  all  were  characterized  by  winds  from  the  north  (Fig.  1)  that  had  not  yet  crossed 
any  hills.  The  northerly  fetch  was  also  the  one  with  the  longest  traverse  over  smooth  snow. 
Although  these  cospectra  have  consistent  behavior  at  high  frequency,  they  too  have  vacillat¬ 
ing  signs  at  low  frequency,  as  in  Figure  7a.  Figures  8b  and  8c  reiterate  the  confused  cospec- 
tral  behavior  at  low  frequency,  with  both  phase  and  coherence  spectra  generally  showing 
poor  correlation  between  u  and  w.  At  higher  frequency,  however,  the  phase  spectra  show 
that  the  turbulence  gets  better  organized;  0u»  is  roughly  180°,  indicating  a  consistent  trans¬ 
fer  of  momentum  to  the  surface  at  these  scales. 

In  summary,  the  topographic  effects  on  the  velocity  field  are  dominant  at  low  frequency. 
Whether  the  flow  is  down  from  the  hills  or  up  over  them,  the  topography  disrupts  the  mo¬ 
mentum  exchange  at  low  frequency,  with  some  scales  transferring  momentum  downward 
while  others  transfer  it  upward.  At  inertial  subrange  frequencies  the  momentum  exchange  is 
still  confused  for  downwind  measurements  within  at  least  6h  of  the  hills,  where  h  is  the 
height  of  the  hills.  For  measurements  upwind  a  distance  of  6h  from  the  hills,  the  inertial  sub¬ 
range  behaves  more  normally,  with  momentum  transfer  toward  the  surface  at  these  scales. 

w-t  and  w-q  cospectra 

Despite  the  topographic  perturbations  the  w-r  and  h-</ cospectra  (Fig.  9a,  lOa,  11a,  12a) 
have  much  more  consistent  behavior  throughout  their  entire  frequency  range  than  the  u-w 
cospectra.  This  difference  may  indicate  that  sampling  times  were  too  short.  Haugen  et  al. 
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c.  Coherence  spectra. 

Figure  8  (cont ’d).  Plots  of  u-w  cospectra,  phase  spectra  and  coherence 
spectra  for  the  runs  in  Table  I  in  which  Co^w  was  negative  in  the  iner¬ 
tial  subrange. 


(1971)  and  Wyngaard  (1973)  recommended  sampling  times  for  stress  of  at  least  one  hour  in 
the  surface  layer  but  found  shorter  times  were  generally  adequate  for  heat  flux.  Because  of 
sporadic  glitches  in  my  recording  system,  I  seldom  had  continuous  runs  longer  than  25-30 
minutes  (Table  1) — about  right  for  sensible  and  latent  heat  fluxes  but,  presumably,  too  short 
for  stress  (Wyngaard  1973). 

Figures  9a,  10a,  1  la  and  12a  contain  cospectra  plotted  nondimensionally  with  typical  sur¬ 
face  layer  scaling  (Kaimal  et  al.  1972).  The  m„  t,  and  q,  values  used  in  the  scaling  are  from 
the  inertial-dissipation  flux  estimates.  I  separated  the  w-t  and  w-q  spectra  into  two  groups 
each,  depending  on  whether  the  cospectra  were  negative  or  positive  in  the  inertial-convective 
subrange.  Because  Figures  9a,  10a,  1  la  and  12a  plot  the  spectra  logarithmically,  they  cannot 
represent  spectral  estimates  with  a  sign  different  from  that  in  the  inertial-convective  sub¬ 
range;  therefore,  some  points  in  each  cospectrum  may  have  been  omitted.  The  dashed  line  in 
each  figure  is  the  Kansas  neutral-stability  w-t  cospectrum  (Kaimal  et  al.  1972): 
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The  majority  of  w-t  cospectra  were  negative  in  the  inertial-convective  subrange  (Fig.  9). 
Despite  the  uncertainties  in  the  inertial-dissipation  technique  over  complex  terrain,  the  de¬ 
rived  u.  and  /,  values  collapse  the  cospectra  well  at  high  frequency  (Fig.  9a):  an  inertial- 
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convective  subrange,  where  the  cospectra  go  as (Wyngaard  and  Cote  1972),  is  fairly  ob¬ 
vious.  Although  the  cospectra  collapse  well,  they  consistently  fall  below  eq  26  at  its  spectral 
peak  and  beyond.  This  suggests  that  the  excess  spectral  energy  at  low  frequency  in  the  u  and  / 
spectra  and  now  in  the  w-t  cospectra  was  indeed  cascading  through  the  inertial  subrange, 
producing  higher  dissipation  rates  than  would  be  found  in  similar  conditions  over  horizon¬ 
tally  homogeneous  terrain.  Consequently  the  uj,  and  u,q,  products  must  be  higher  (see  eq 
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a.  w-t  cospectra. 


b.  Phase  spectra. 

Figure  9.  Plots  of  w-t  cospectra,  phase  spectra  and  coherence  spectra  for 
16  runs  in  which  Co^,,  was  predominantly  negative  in  the  inertial-convec¬ 
tive  subrange.  The  dashed  line  is  the  Kansas  neutral-stability  w-t  cospectrum 
(eq  26). 
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c.  Coherence  spectra. 

Figure  9  (cont‘d).  Plots  of  w-t  cospectra,  phase  spectra  and  coherence 
spectra  for  16  runs  in  which  COnif  was  predominantly  negative  in  the  iner¬ 
tial-convective  subrange.  The  dashed  line  is  the  Kansas  neutral-stability  w-t 
cospectrum  (eq  26). 


10-12).  This  result  is  also  evident  in  Figure  9a.  Because  the  integral  of  COwt  yields  mT  = 
integrating  the  measured  w-t  cospectra — with  their  excess  energy  at  low  frequency — 
will  produce  larger  ti7  values  than  integrating  the  Kansas  w-t  cospectrum  (the  dashed  line  in 
Fig.  9a). 

The  phase  and  coherence  spectra  (Fig.  9b  and  c)  emphasize  the  more  consistent  behavior 
of  the  w-t  cospectra.  For  nondimensional  frequencies  above  0. 1 ,  the  phase  spectra  cluster 
around  180°,  indicating  that  the  w  and  t  fluctuations  had  good  negative  correlation.  At  low¬ 
er  frequencies  the  phase  spectra  are  more  scattered,  but  the  bulk  of  the  points  are  between 
90°  and  270°.  That  is,  the  correlation  generally  remains  negative  at  low  frequency:  the  spec¬ 
tral  estimates  outside  this  range  indicate  that  some  cospectral  estimates  have  gone  positive. 
Compare  the  u-w  phase  spectra  (Fig.  7b).  At  high  frequency,  ©u*  is  typically  between  0  °  and 
90  °  (Couw  positive),  while  at  low  frequency  the  majority  of  the  spectral  estimates  are  in  the 
180-270°  quadrant  (COy*  negative).  The  w-t  coherence  spectra  (Fig.  9c)  show  moderate  cor¬ 
relation  at  low  frequency.  With  the  high  energy  levels  in  the  w  and  t  spectra  here,  this  results 
in  the  high  cospectral  energy  evident  in  Figure  9a. 

The  w-t  cospectra  that  are  predominantly  positive  in  the  inertial-convective  subrange  are 
more  erratic  (Fig.  10).  The  collapse  of  the  cospectra  when  scaled  with  u,t,  is  not  very  success¬ 
ful,  even  at  high  frequency  (Fig.  10a).  The  phase  spectra  (Fig.  10b)  show  that  the  sign  usually 
changes  between  low  (negative)  and  high  (positive)  frequencies.  The  coherence  spectra  (Fig. 
10c)  are  not  dramatically  different  from  those  in  Figure  9c  except  in  the  extended  portion  of 
the  spectrum,  where  a  higher  percentage  of  spectral  estimates  have  Chy,^  values  much  larger 
than  I  (and,  hence,  off  the  plot)  than  for  the  negative  cospectra.  The  probable  cause  of  this 


erratic  behavior  is  the  low  energy  in  the  positive  w-t  cospectra.  Table  2  shows  that  the  runs 
with  COwt  positive  in  the  inertial-convective  subrange  (positive  H^)  had  some  of  the  smallest 
absolute  values  of 

An  overwhelming  majority  of  the  w-q  cospectra  are  positive  in  the  inertial-convective 
subrange  (Fig.  1  la).  Again  the  surface-layer  nondimensionalization  collapses  these  cospectra 
well  at  peak  frequencies  and  beyond.  An  inertial-convective  subrange,  where  Co^q  goes  as 
is  evident.  But,  as  with  the  w-t  cospectra,  at  high  frequency  the  nondimensional  co¬ 
spectra  fall  below  the  neutral-stability  Kansas  w-t  cospectrum  (eq  26).  I  attribute  this  again 
to  higher  values  of  u,q,  at  this  topographically  complex  site  than  would  derive  over  horizon¬ 
tally  homogeneous  terrain.  At  low  frequency  the  w-q  cospectral  estimates  are  not  as  consis¬ 
tently  above  the  Kansas  cospectrum  as  were  the  w-t  cospectral  estimates  in  Figure  9a.  This 
could  be  a  consequence  of  poorer  coherence  between  the  w  and  q  fluctuations  (Fig.  1  Ic)  or, 
more  probably,  the  slightly  lower  energy  in  4>q  than  in  at  low  frequency. 
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a.  w-t  cospectra. 
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b.  Phase  spectra. 

Figure  W.  Plots  of  w-t  cospectra,  phase  spectra  and  coherence  spectra  for 
five  runs  in  which  Co„i  was  positive  in  the  inertial-convective  subrange.  The 
dashed  tine  is  the  Kansas  neutral-stability  w-t  cospectrum  (eq  26). 
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c.  Coherence  spectra. 

Figure  10  (cont’d).  Plots  of  w-t  cospectra,  phase  spectra  and  coherence 
spectra  for  five  runs  in  which  Co^,,  was  positive  in  the  inertial-convective 
subrange.  The  dashed  line  is  the  Kansas  neutral-stability  w-t  cospectrum  (eq  26). 


a.  w-q  cospectra. 


Figure  11.  Plots  of  w-q  cospectra,  phase  spectra  and  coherence  spectra  for 
22  runs  in  which  Co„g  was  positive  in  the  inertial-convective  subrange.  The 
dashed  tine  is  the  Kansas  neutral-stability  w-t  cospectrum  (eq  26). 
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6.  Phase  spectra. 
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c.  Coherence  spectra. 


Figure  1 1  (cont’d). 


The  cospectrum  and  phase  and  coherence  spectra  for  the  one  run  (11.2)  for  which  Co^^q 
was  negative  in  the  inertial-convective  subrange  are  shown  in  Figure  12.  Again  COv»q  falls  be¬ 
low  the  Kansas  w-t  cospectrum  at  its  peak  and  beyond.  The  phase  spectrum  indicates  consis¬ 
tent  negative  correlation  at  high  frequency  but  sign  reversals  at  low  frequency. 

The  vacillating  sign  at  low  frequency  and  the  frequent  sign  difference  between  low  and 
high  frequencies  in  the  w-t  and  w-q  cospectra  suggest  little  coherence  between  transfer  pro¬ 
cesses  at  large  and  small  scales  at  our  site.  The  small  scales  seem  to  be  near  equilibrium  with 
the  snow  surface  because  inertial-dissipation  estimates— based  on  the  high-frequency  behav¬ 
ior  of  4>u,  <t>,  and  <I>q  and  on  the  signs  of  COw,  and  Co^q  in  the  inertial-convective  subrange— 
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yield  reasonable  values  of  and  and  values  of  u„  t,  and  q,  that  are  effective  in  collaps¬ 
ing  Co*t  and  Co^^q  to  consistent  form.  These  values  are  also  consistent  with  upper-air 
soundings  made  during  the  experiment.  Although  soundings  were  not  always  made  twice  a 
day  during  the  experiment,  the  available  soundings  consistently  showed  a  shallow  inversion 
layer,  with  a  well-mixed  layer  above.  Therefore,  the  predominantly  negative  values  that 
the  inertial-dissipation  technique  yielded  (Table  2)  are  consistent  with  atmospheric  boundary 
layer  profiles. 

t-q  cospectra 

Because  the  w-t  cospectra  were  generally  negative  in  the  inertial-convective  subrange 
while  the  w-q  cospectra  were  virtually  always  positive,  the  temperature  and  humidity  fluctu¬ 
ations  were  often  negatively  correlated.  The  t~q  cospectra  reiterate  this  conclusion  and  point 
out  other  interesting  aspects  of  temperature  and  humidity  transfer  over  complex  terrain. 
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a,  w-q  cospectrum. 


b.  Phase  spectrum. 

Figure  12.  Plots  of  the  w-q  cospectrum,  phase  spectrum  and  co¬ 
herence  spectrum  for  the  one  run  (1 1.2)  in  which  COy„q  was  nega 
tive  in  the  inertial-convective  subrange.  The  dashed  line  is  the  Kan¬ 
sas  neutral-stability  w-t  cospectrum  (eq  26). 
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c.  Coherence  spectrum. 
Figure  12  (cont’d). 


As  with  the  w-t  and  w-q  cospectra,  I  have  grouped  the  t-q  cospectra  according  to  whether 
they  are  negative  or  positive  in  the  inertial-convective  subrange.  Figure  13  shows  cospectra 
and  phase  and  coherence  spectra  for  t-q  cospectra  that  are  negative  in  the  inertial-convective 
subrange;  Figure  14  shows  the  same  plots  for  t-q  cospectra  that  are  positive.  The  nondimen- 
sionalization  in  Figures  13a  and  I4a  again  reflects  surface-layer  scaling  based  on  the  predic¬ 
tion  that  Cotq  should  behave  in  an  inertial-convective  subrange  just  like  temperature  and  hu¬ 
midity  spectra  (Friehe  et  al.  1975,  Wyngaard  et  al.  1978): 

Co^^ik)  =  (27) 

Here  is  the  dissipation  rate  of  the  t-q  covariance;  and  I  assume  0,q  =  0.40,  half  the  value 
used  for  the  temperature  and  humidity  spectra  (Antonia  et  al.  1978,  Hill  1978).  Since  Co,q  is 
assumed  to  behave  like  temperature  and  humidity  spectra,  the  dashed  line  in  Figures  13a  and 
14a  is  the  Kansas  neutral-stability  temperature  spectrum  (eq  6). 

The  surface-layer  nondimensionalization  works  well  for  both  negative  and  positive  t-q  co¬ 
spectra.  The  cospectra,  especially  the  negative  t-q  cospectra,  collapse  well  in  the  inertial- 
convective  subrange  (Fig.  13a  and  14a).  Both  sets  of  cospectra,  however,  generally  have 
more  energy  than  the  Kansas  temperature  spectrum  for  frequencies  at  its  peak  and  lower.  Be¬ 
cause  of  this  excess  energy  at  low  frequency,  the  t-q  cospectra  have /  regions  that  begin  at 
a  frequency  roughly  an  order  of  magnitude  lower  than  in  the  Kansas  spectrum.  Remember, 
we  observed  these  same  features  in  the  temperature  and  humidity  spectra.  Thus,  although 
my  scalar  spectra  and  cospectra  do  not  look  exactly  like  the  Kansas  temperature  spectrum, 
they  look  very  much  like  each  other. 
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b.  Phase  spectra. 

Figure  13.  Plots  of  t-q  cospectra,  phase  spectra  and  coherence 
spectra  for  12  runs  in  which  Co,q  was  predominantly  negative  in  the 
inertial-convective  subrange.  The  dashed  line  is  the  Kansas  neutral- 
stability  temperature  spectrum  (eg  6). 
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c.  Coherence  spectra. 
Figure  13  (cont’d). 


a.  t-q  cospectra. 


Figure  14.  Plots  of  t-q  cospectra,  phase  spectra  and  coherence  spec¬ 
tra  for  five  runs  in  which  Co,^  was  positive  in  the  inertial-convective 
subrange.  The  dashed  line  is  the  Kansas  neuiral-stahiliiy  temperature  spec¬ 
trum  (eq  6). 
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b.  Phase  spectra. 


c.  Coherence  spectra. 

Figure  14  (cont’d).  Plots  of  t-q  cospectra,  phase  spectra  and 
coherence  spectra  for  five  runs  in  which  Co,q  was  positive  in  the  iner¬ 
tial-convective  subrange.  The  dashed  line  is  the  Kansas  neutral-stability 
temperature  spectrum  (eq  6). 


The  one  anomalous  feature  of  the  t-q  cospectra  is  that  the  inertial-convective  subrange, 
the  f  region,  does  not  extend  to  frequencies  as  high  as  it  does  in  the  t  and  q  spectra,  al¬ 
though  it  theoretically  should.  The  reason  is  the  separation  between  the  t  and  q  sensors,  d  = 
4  cm.  Kaimal  (1975)  explained  that  as  the  turbulence  frequency  gets  higher  and  the  eddy  size 
becomes  comparable  to  the  separation  between  sensors,  the  turbulence  at  the  two  sensors  be¬ 
comes  increasingly  uncorrelated.  The  nondimensional  frequency /z/f/ at  which  this  occurs  is 
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roughly  z/6ird.  For  my  measurements  this  value  is  about  3  [or  log(z/67r^/)],  which  corre¬ 
sponds  well  with  the  frequency  in  Figures  13a  and  14a,  where  the  cospectra  turn  down  from 
the  -V}  slope.  The  phase  spectra  (Fig.  13b  and  14b)  also  start  losing  their  inertial-subrange 
consistency  near  this  frequency  (Priestley  and  Hill  1985).  Recognizing  this  effect  in  the  t-q 
cospectra,  we  can  also  go  back  and  see  it  in  the  more  scattered  w-t  and  w-q  cospectra.  The 
hot-film  sensors  that  yielded  w  were  4  cm  from  the  Lyman-alpha  hygrometer  and  2.5  cm 
from  the  PRT.  Hence,  the  change  in  slope  for  Co*,  and  Co*q  in  Figures  9a  and  1  la  occurs  at 
roughly  the  same  frequency  as  for  Co,q. 

Because  few  measurements  of  t-q  cospectra  have  been  reported,  the  validity  of  the  inertial- 
convective  subrange  form  of  the  cospectrum  (eq  27)  is  still  tentative.  Friehe  et  al.  (1975)  pre¬ 
sented  one  negative  cospectrum  obtained  over  the  Salton  Sea  and  one,  likely  positive,  co¬ 
spectrum  measured  over  the  ocean.  Both  had  the  predicted  inertial-convective  subrange  be¬ 
havior.  Wyngaard  et  al.  (1978)  showed  one  positive  t-q  cospectrum  that  behaved  as  in 
the  inertial-convective  subrange  and  said  it  was  typical  of  the  40  cospectra  that  they  had  ob¬ 
tained  over  the  East  China  Sea.  They  also  stated  that  unpublished  t-q  cospectra  computed 
from  data  collected  by  Champagne  et  al.  (1977)  over  Minnesota  farmland  also  showed  iner¬ 
tial-convective  subranges.  These  were  likely  also  positive  cospectra.  Antonia  and  Chambers 
(1980)  also  showed  one  positive  cospectrum  that  went  as at  high  frequency  and  was  rep¬ 
resentative  of  the  four  they  had  measured  over  Bass  Strait.  Priestley  and  Hill  (1985)  meas¬ 
ured  six  positive  and  six  negative  t-q  cospectra  at  BAO  that  behaved  as  McBean  and 
Elliott  (1981)  and  Kohsiek  (1982),  however,  did  not  find  the  predicted  inertial-convective 
subrange  behavior  in  the  t-q  cospectra  that  they  obtained  over  dry  land.  Figure  13a  is  thus 
important  confirmation  of  the  existence  of  an  inertial-convective  subrange  in  negative 
t-q  cospectra. 

As  with  essentially  all  the  cospectra  that  I  have  presented,  the  t-q  cospectra  tend  to  have 
different  signs  at  low  and  high  frequencies.  The  phase  spectra  (Fig.  13b  and  14b)  show  this 
clearly.  For  the  runs  with  Co,q  positive  in  the  inertial-convective  subrange,  this  tendency  is 
not  as  pronounced  as  when  Cotq  is  negative.  Nevertheless  the  phase  spectra  suggest  that  the 
scalar  content  of  large  and  small  eddies  was  often  different.  For  example,  for  Co,q  negative 
in  the  inertial-convective  subrange  (Fig.  13),  at  high  frequency  warm  eddies  were  generally 
drier  than  ambient  conditions,  while  cool  eddies  were  more  moist.  At  low  frequency,  on  the 
other  hand,  warm  eddies  were  often  more  moist  than  ambient  conditions,  while  cool  eddies 
were  drier. 

The  t-q  coherence  spectra  (Fig.  13c  and  14c)  are  different  from  any  of  the  other  coherence 
spectra.  First,  the  coherence  is  higher  at  high  frequency  than  for  u-w,  w-i  or  w-q.  This  is  a 
consequence  not  only  of  the  better  correlation  (both  positive  and  negative)  between  tempera¬ 
ture  and  humidity  fluctuations  at  high  frequency  but  also  of  the  less-rapid  roll-off  of  Co,q 
than  of  COu*,  Co*,  and  Co*q  (Fairall  et  al.  1980).  The  one  run  that  had  C/j,q  values  conspic¬ 
uously  above  0.90  in  Figure  14c  was  1 1 .2 — the  only  run  during  which  and  f/,  were  both 
negative  (Table  2).  This  seems  to  suggest  that  t  and  q  have  a  better  positive  correlation  at 
high  frequency  in  stable  conditions  than  in  unstable  conditions.  Although  few  corroborating 
coherence  spectra  have  been  reported,  Priestley  and  Hill  (1985)  did  point  out  that  in  their 
data — for  which  was  always  positive — the  t-q  coherence  was  usually  higher  in  stable  than 
in  unstable  conditions. 

A  second  important  feature  of  the  coherence  spectra  is  the  dichotomy  indicated  in  Figure 
13c.  The  coherence  spectra  are  generally  high  at  low  frequency,  have  local  minima  near  fz/U 
=  0.02,  then  peak  again  in  the  inertial-convective  subrange.  The  conclusion  must  again  be 
that  large-scale  and  small-scale  transport  processes  are  dissimilar.  The  phase  spectra  in  Fig¬ 
ure  13b  corroborate  this  conclusion;  they  are  consistently  near  180°  above  fz/V  =  0.02  and 
scattered  at  lower  frequencies.  The  runs  that  yielded  Figure  13  were  all  characterized  by 


From  eq  28  the  refractive  index  spectrum  ^n{k)  at  wavenumber  k  is 


<l>n(A:)  =  /!  =  <*>, (A:)  +  B^<^^(k)  +  2ABCo^^(k).  (29) 

Therefore,  by  making  simultaneous  measurements  of  the  temperature  and  humidity  spectra 
and  the  t-q  cospectrum,  I  can  compute  the  refractive  index  spectrum  for  an  arbitrary  electro¬ 
magnetic  wavelength.  Since  4>t(A:),  'l>q(A:)  and  Co,q(A:)  are  all  proportional  to  -  ’  in  an 
inertial-convective  subrange,  4>n(^)  should  be  also.  Specifically, 

<l>n(A:)  = (30) 

where  is  the  dissipation  rate  of  half  the  refractive  index  variance,  and  (i  is  the  same  Kol¬ 
mogorov  constant  I  have  used  for  the  other  scalar  spectra  (Hill  1978). 

In  studies  of  electromagnetic  wave  propagation  in  the  atmosphere,  another  representation 
for  spectra  in  the  inertial-convective  subrange  arises.  This  involves  the  structure  parameters 
Q,  q,  C,q  and  q: 


4>t(A:)  =  0.249q  k'^'^ 

(31) 

<J>q(A:)  =  0.249qA:-5^^ 

(32) 

Co, q(A:)  =  0.249Ctq  k'^'^ 

(33) 

4>n(A:)  =  0.249C^A:-5''\ 

(34) 

Clearly,  from  eq  29  the  refractive  index  structure  parameter  is  related  to  the  structure  param¬ 
eters  for  t,  q  and  t-q  by 

C^n  =  +  B^C^  +  2ABQ^.  (35) 

Scintillometers,  for  example,  measure  directly.  Thus,  two  ways  exist  to  investigate  the  re¬ 
fractive  index  spectrum — with  electro-optical  instruments  or  with  small,  fast-responding  me¬ 
teorological  instruments.  Portman  et  al.  (1962)  and  Ryznar  (1963)  did  some  of  the  early 
work  on  electro-optical  propagation  over  snow  and  its  relation  to  atmospheric  turbulence. 

Figures  15  and  16  show  0.55-^im  ('t’sj)  and  millimeter  refractive  index  spectra  com¬ 

puted  using  eq  29  from  the  temperature  and  humidity  spectra  and  the  l-q  cospectra  that  I 
have  been  discussing.  The  spectra  are  again  made  nondimensional  with  surface-layer  scaling, 
and  the  dashed  line  is  the  modified  Kansas  temperature  spectrum  (eq  6).  These,  1  believe,  are 
the  first  refractive  index  spectra  for  either  0.55-tim  or  millimeter  wavelength  light  ever  ob¬ 
tained  over  snow.  And  only  a  few  visible  (Friehe  et  al.  1975,  Antonia  et  al.  1978)  and  milli¬ 
meter  (Priestley  and  Hill  1985)  wavelength  refractive  index  spectra  have  been  published  for 
other  surfaces. 

The  two  sets  of  refractive  index  spectra  look  like  each  other  and  also  like  the  scalar  spectra 
and  cospectra  1  presented  earlier.  That  is,  both  sets  of  spectra  have  more  energy  at  low  fre¬ 
quency  than  does  the  Kansas  neutral-stability  temperature  spectrum — again  because  of  topo¬ 
graphic  effects.  Both  sets,  however,  join  the  Kansas  spectrum  at  higher  frequency  but  evi¬ 
dently  begin  their  inertial-convective  subranges  at  lower  frequency  than  in  the  Kansas  spec¬ 
trum. 
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Table  4.  Values  of  the  spectrally  derived  re¬ 
fractive  index  structure  parameters  for 
0.55-^  and  millimeter  wavelengths  and 
the  relative  contributions  of  temperature 
and  humidity  fluctuations  and  the  t-q  co- 
variance  to  the  millimeter  structure  param- 


eter. 

10"  C;(m-^’7 

Contributions  to  mm  Cj^ 

Run 

0.55  itm 

mm 

t 

9 

t-q 

l.l 

3.57 

9.57 

0.38 

0.29 

0.33 

1.2 

2.19 

9.38 

0.24 

0.41 

0.35 

2.1 

2.44 

10.1 

0.24 

0.39 

0  36 

2.2 

2.96 

14.4 

0.21 

0  39 

0.40 

3.1 

0.18 

0.78 

0.23 

0.73 

0.04 

4.1 

0.09 

3.76 

0.02 

1.07 

-0.09 

4.2 

0.14 

5.37 

0.03 

0.86 

0.11 

5.1 

0.19 

15.1 

0.01 

0.88 

0.11 

5.2 

0.14 

15.2 

0.01 

0.90 

0.09 

6.1 

2.29 

6.2 

3.14 

6.3 

2.46 

7.1 

1.97 

7.2 

2.66 

9.1 

0.26 

9.07 

0.03 

0.73 

0.24 

9.2 

0.39 

9.58 

0.04 

0.71 

0.25 

10.1 

0.47 

11.7 

0.04 

0.70 

0.26 

10.2 

0.25 

6.37 

0.04 

0.69 

0.27 

10.3 

1.44 

12.2 

0.12 

0.70 

0.17 

11.1 

3.51 

3.55 

1.01 

0.09 

-0.10 

11.2 

6.36 

2.24 

2.88 

0.80 

-2.67 

12.1 

0.17 

7.43 

0.02 

1.09 

-o.n 

13.1 

2.06 

2.20 

0.94 

0.13 

-0.07 

13.2 

4.72 

5.55 

0.85 

0.07 

0.07 

14.1 

2.61 

2.37 

1.10 

0.45 

-0.55 

With  the  presence  of  a  conspicuous  inertial-convective  subrange  in  each  refractive  index 
spectrum,  it  is  possible  to  find  Q  for  each.  Table  4  lists  these  values.  The  table  also  lists  the 
contributions  to  the  millimeter  value  from  the  temperature  and  humidity  spectra  and  the 
t-g  cospectrum  (see  eq  35).  Despite  the  low  temperatures,  humidity  fluctuations  generally 
make  the  largest  contribution  to  the  millimeter  Cp  value.  The  t-q  covariance  is  usually  next 
in  importance  and  can  make  a  negative  contribution,  thereby  offsetting  relatively  large  t  or  q 
contributions.  The  temperature  fluctuations,  though  generally  making  the  smallest  contribu¬ 
tion,  are  seldom  negligible.  Consequently  all  three  terms  in  eq  29  or  35  are  necessary  for  ac¬ 
curately  measuring  <!>„  or  in  the  millimeter  region  over  snow. 

The  temperature  spectrum  always  contributed  at  least  96%  to  the  0.55-^im  value,  so  1 
have  not  provided  a  similar  breakdown  for  this  wavelength  in  Table  4.  McBean  and  Elliott 
(1981)  also  found  that  temperature  fluctuations  far  outweighed  other  contributions  to  in 
the  visible  region  over  land.  Over  the  ocean,  on  the  other  hand,  Antonia  et  al.  (1978)  report¬ 
ed  that  the  t-q  contribution  to  Cp  could  be  5-10%  even  though  the  q  contribution  was  less 
than  1%. 

During  the  experiment,  another  group  measured  Cp  for  visible  wavelengths  with  an  optical 
scintillometer  (Wang  et  al.  1978)  placed  near  my  instrument  mast.  The  optical  path  was  300  m 
long  and  roughly  2  m  above  the  surface.  Figure  17  compares  my  spectrally  derived  0.55-/(m 
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Figure  17.  Comparison  of  simultaneous 
scintillometer  and  spectrally  derived  values 
of  C^nfor  visible  light. 

Q  values  with  these  simultaneous  scintillometer  measurements.  The  two  methods  agree  rea¬ 
sonably  well  at  high  spectral  but  not  as  well  at  low  spectral  Q.  Seeing  this  and  then  look¬ 
ing  more  closely  at  how  the  scintillometer  was  set  up,  1  found  that  it  was  operating  too  near 
its  noise  level  (Ochs  et  al.  1977)  and,  thus,  was  incapable  of  resolving  the  low  Cl^  values  en¬ 
countered.  Hence,  this  first  experimental  attempt  to  verify  eq  35  over  snow  was  not  entirely 
successful.  Kohsiek  (1985),  however,  did  find  good  agreement  between  infrared  Q  values 
obtained  with  a  COj  laser  and  with  fast-responding  meteorological  instruments  over  varied 
terrain  in  the  Netherlands.  And  Priestley  and  Hill  (1985)  found  reasonable  agreement  be¬ 
tween  radio-wavelength  Q  values  obtained  at  BAO  with  a  radio  refractometer  and  with  me¬ 
teorological  instruments.  Evidence  is,  therefore,  mounting  that  the  theoretical  arguments 
that  lead  to  eq  29  and  35  are  basically  sound. 


CONCLUSIONS 


1  have  described  an  extensive  set  of  turbulence  spectra  and  cospectra  collected  in  near¬ 
neutral  conditions  over  snow.  All  of  the  spectra — u,  w,  t  and  q —  and  cospectra — u-w,  w-t, 
w-q  and  t-q — have  high  energy  at  low  frequency  because  of  flow  disturbances  associated 
with  the  forested  hills  bordering  the  site.  With  the  exception  of  Co^^,,  all,  nevertheless,  have 
the  expected  inertial  or  inertial-convective  subrange  behavior.  The  erratic  behavior  of  Coyw 
is,  thus,  primary  evidence  of  flow  disturbance. 

The  phase  and  coherence  spectra,  especially  Chtq,  suggest  a  difference  in  transfer  process¬ 
es  at  large  and  small  scales  at  this  topographically  complex  site.  With  the  atmospheric  boun¬ 
dary  layer  generally  having  slightly  stable  stratification,  1  found  evidence  of  internal  gravity 
waves  when  the  surface  layer  was  also  stably  stratified.  These  low-frequency  waves  and  the 
low-frequency  turbulence  induced  by  the  topography  generally  combined  to  produce  con¬ 
fused  cospectra  with  vacillating  sign  at  low  frequency.  At  higher  frequency  the  spectra  and 
cospectra  were  better  behaved,  thus  indicating  that  here  the  fluctuations  resulted  strictly 
from  turbulent  processes  that  were  nearly  in  equilibrium  with  the  smooth,  flat  snow  surface 
surrounding  my  instrument  mast. 
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With  simultaneous  measurements  of  temperature  and  humidity  spectra  and  temperature- 
humidity  cospectra,  I  could  compute  refractive  index  spectra  for  light  of  visible  and  milli¬ 
meter  wavelengths.  These  refractive  index  spectra  looked  like  the  scalar  spectra  that  yielded 
them;  they  had  significant  energy  at  low  frequency  and  an  inertial -convective  subrange  at 
high  frequency.  Over  snow  at  the  air  temperatures  encountered,  the  spectrum  for  visible  light 
is  determined  almost  entirely  by  temperature  fluctuations.  For  the  millimeter  spectrum,  on 
the  other  hand,  humidity  fluctuations  are  generally  most  important,  but  temperature  fluctu¬ 
ations  and  the  t-q  covariance  cannot  be  ignored.  These  spectra  yielded  refractive  index  struc¬ 
ture  parameters  Q  that  ranged  from  9  x  10  ”  to  6  x  lO"”  m"^’  for  visible  light  and  from  8  x 
lO'"  to  1.5  X  lO'”  m'^‘  for  millimeter  wavelengths. 
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